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AN ELECTROSTATIC INTERPRETATION OF THE ZEROS OF 
PARAORTHOGONAL POLYNOMIALS ON THE UNIT CIRCLE 

BRIAN SIMANEK 


Abstract. We show that if /i is a probability measure with infinite support on the unit 
circle having no singular component and a differentiable weight, then the corresponding 
paraorthogonal polynomial <& n (z\fi) solves an explicit second order linear differential equa¬ 
tion. We also show that if r 7 ^ /3, then the pair (& n (z\ /3), <£> n (,s; r)) solves an explicit first 
order linear system of differential equations. One can use these differential equations to 
deduce that the zeros of every paraorthogonal polynomial mark the locations of a set of 
particles that are in electrostatic equilibrium with respect to a particular external field. 
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1. Introduction 

Given a positive probability measure /i with compact and infinite support in the com¬ 
plex plane C, one can constuct the corresponding sequence of orthonormal polynomials 
{< p n (z ; /i)} n >o where <p n (z‘, /i) is a polynomial of degree exactly n having positive leading co¬ 
efficient, which we denote by K. n (/i). If we divide ip n (z; //) by re n (//), then we obtain a monic 
degree n polynomial, which we denote by <3> n (s; /i). We will often suppress the /j dependence 
in our notation if there is no possibility for confusion. Our study of these objects is motivated 
by the fact that orthogonal polynomials have proven to be valuable tools in the study of 
physical models. 

The most well-studied collections of orthogonal polynomials come from measures of or¬ 
thogonality supported on the real line and are known as orthogonal polynomials on the real 
line (OPRL). Such polynomials have a variety of applications in spectral theory (see |3T]), 
potential theory (see Pi). and the theory of special functions (see na). One of the key 
features of OPRL is the three-term recurrence formula satisfied by the orthonormal polyno¬ 
mials. Many well-known families of OPRL (including the Hermite polynomials, the Laguerre 
polynomials, and the Jacobi polynomials) are polynomial solutions to particular families of 
linear second order differential equations with polynomial coefficients. This property makes 
these particular families of polynomials especially applicable as we will discuss later in more 
detail. 

A second collection of well-studied orthogonal polynomials comes from measures of or¬ 
thogonality supported on the unit circle and are known as orthogonal polynomials on the 

unit circle (OPUC). Just as OPRL has a close connection with self-adjoint operators, OPUC 
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has a close connection with unitary operators (see [3T]). In analogy with the three-term 
recurrence satisfied by the OPRL, the monic OPUC satisfy the Szego recursion: 

<P n+ i(z; p) = z$ n (z] p) - a n $*(z; p), 

where a n e D := {z : \z\ < 1} and <I>*(z) = z n Q n (l/z). Therefore, to each infinitely sup¬ 
ported probability measure p on the unit circle, we can associate the sequence {a n }n>o of 
Verblunsky coefficients. A theorem of Verblunsky states that such a sequence also deter¬ 
mines an infinitely supported probability measure (see [281 Chapter 1]). The utility and 
applicability of the Hermite, Laguerre, and Jacobi polynomials inspired interest in families 
of OPUC that are solutions to linear second order differential equations. To this end, Ismail 
and Witte proved the following result in [ 16] : 


Theorem 1.1 (Ismail & Witte, 2001). Letw(z) = e v(z ' > be differentiable in a neighborhood 
of the unit circle, have moments of all integral orders, and assume that the integrals 



v'(z) — v'(() 
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exist for all integers n. Then the corresponding orthonormal polynomials satisfy the differ¬ 
ential relation 

<M) = A-ifC - MO) [ —t-^MCMCMCM 

K n •'ICI =1 Z — C 

+ iifn(z) f oil) — ^Q-(p n (()<p n (()w(()d(. 

J|C|=i z ~ C 


From this result, one can derive a second order differential equation to which <p n (z) is a 
solution; we refer the reader to [IB] for details. The main purpose of our investigation is to 
derive a comparable result for a related class of polynomials called paraorthogonal polynomials 
on the unit circle (POPUC). 

Given an infinitely supported probability measure p on the unit circle and a complex 
number fd of modulus 1, we define the paraorthogonal polynomial Q n (z] fd] p) as the monic 
degree n polynomial given by 


MuM) := £$ n -i(z;/x) - P$* n -i(z]p)- 


( 1 ) 


Paraorthogonal polynomials and their zeros have received considerable recent attention from 
the research community (see for example [2j j8[ FT. HE, US, 'SH E2, |35[ 37]). It is well- 
known and easy to show that all of the zeros of $ n (x;; /3; p) are distinct and lie on the unit 
circle. Furthermore, if r ^ fd are distinct complex numbers of modulus 1, then the zeros of 
jd ; p) and & n (z] r; p) strictly interlace on the unit circle in that if x and y are two zeros 
of $ n (^; fd] p) and [x,y\ is the arc of the unit circle that runs from x to y in the counter¬ 
clockwise direction, then [x, ?/]\{x, y} contains a zero of <1> n (z; r; p). Although paraorthogonal 
polynomials are not orthogonal polynomials, they often serve as an appropriate analog of 
OPRL in settings where the real line is replaced by the unit circle (see for example [18,135]). 
One of our main results is the next theorem, which can be thought of as an analog of the 
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Ismail and Witte result that applies to paraorthogonal polynomials. In fact, the theorem 
applies to any degree n monic polynomial satisfying the relation (Jl| for some complex number 
j3. As in [28], we write a measure on the unit circle as a measure in the variable 6 G 
[0, 27t) . For brevity and because there is no possibility for confusion, we will suppress the 
/r-dependence of various quantities in our notation. 


Theorem 1.2. Suppose dy{8) = w{6 )|^ is a probability measure on the unit circle, where 
w is continuous on [0,27r] (mod 2ir) and differentiable on (0, 27t) and let {a n }ff=o be the 
corresponding sequence of Verblunsky coefficients. If (3 E C, then the polynomial y(z ) = 
$> n {z\f$) defined by 0 solves the following differential equation on any domain including 
infinity or zero on which the coefficients are meromorphic: 


0 = y"{z) + 

+ 


1 -n _ h' n {z; (3-(3) 
z h n (z;/3;/3) 


y'(z) 


W[h n (z;f3;f3),h n (z;-P;P)] 1 \\ 

oh i i r r\ ((n + zG n (z))G n (z) + J n {z)[D n (z) na. n _ i)) 

2f3zh n (z-,/3] (3) z 


where 


( 2 ) 

y{z), 


r2 * \Vn-i(e ie )\ 2 w'{8) dd 
(z — e ld ) 27r 

( Pn-i(e id ) 2 w'(6) d6 

o (z — e ld )e ind 2 tt 

^n-1(eg)Mg) dd 

(z — e ie )ed n ~ 2 ) e 2n 


GJz) : = i 


c2tt 


D n (z ) : = -%z 


Jni^z) ■ i 


h n {z ; x; y) : = x(n( 1 - ya n -f) + zG n [z ) + yD n [z )) - z(J n (z) - yG n (z)), 


and W[f,g\ denotes the Wronskian of f and g (that is, W[f,g] = fg' — gf). 


By allowing the parameter f3 to be arbitrarily chosen in C, we can apply our results to 
a wide variety of polynomials including paraorthogonal polynomials and perturbations of 
orthogonal polynomials without a precise understanding of the induced perturbation to the 
measure of orthogonality. 

In our applications of Theorem 1.2, we will focus on the case \f3\ — 1 so that the polynomial 
$ n (^; (3] fj) is a paraorthogonal polynomial. Notice that if we set f3 = a n - 1 , then <3> n (z; (3] y) 
is just the monic orthogonal polynomial < F ri (A; y). In this case, Theorem 1.2 yields a second 
order differential equation solved by $ n (z;//). If we set (3 — 0, then Theorem 1.2 yields 
a second order differential equation solved by z$ n -\(z‘, y), which can be rewritten as a 
differential equation solved by $ n _i (z; y). There is no reason to think that the second order 
differential equations solved by < F ri (^; y) derived by these two methods will be the same as 
each other or the same as the one produced by the results in [T(J if those results apply. In 
fact, we will see by example in Section |4.4 that these methods may yield different differential 
equations. It is easy to see that the polynomial y) satisfies many differential equations 
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if one allows the complexity of the coefficients to grow with n. Indeed, if 

K& V) + P { Z Wn( Z ', I*) + Q(z)$n(z] fi) = 0, 
and R(z) is any entire meromorphic function, then 

K( z ; m) + K( z ; v) ( p i z ) + R(z)$ n (z; p)) + $ n (z; p) ( Q{z ) - R(z)& n (z; fi)) = 0. 


In the next section we will discuss applications of Theorem L2 to electrostatics. In that 
context, the coefficients in the second order differential equation determine the location of 
charges that generate an electric field and the zeros of <I> n (z; /3; fi) mark the location of charges 
in equilibrium with respect to that field. It is not surprising that a particular configuration 
of charges can be in equilibrium with respect to different electric fields, so from a physical 
perspective, the lack of uniqueness of the differential equation is also expected. 


An intermediate step in the proof of Theorem |1.2| is of interest in its own right and can 
be stated as follows: 


Theorem 1.3. Suppose p is as in Theorem \1.2\ and t,/ 3 are distinct complex numbers. 
The polynomials u(z) := <S> n (z\ fd, p) and v(z ) := & n (z] r, fi) solve the following system of 
differential equations on any domain containing infinity or zero on which the coefficients are 
meromorphic: 


,/ n / w X 0 ;/?;/ 3 )\ / x fhn(z;T‘,P) 

U (z) = v[z) - ~r= -— — U{z) 1 


v'(z) = v(z) 


z (/3 - T) 

h n (z;/3;r) 


u(z ) 


z{P ~ T ) 

h n (z ; r; r) 


( 3 ) 

( 4 ) 


z{P - r ) ) V z{P -t) 

The main restriction in the applicability of Theorems L2 and IH is the requirement that 


the measure be given by a continuous and differentiable weight function. However, for any 
fixed neN and /3 G <90, the map from measures to degree n paraorthogonal polynomials 
that map 0 to ft is far from injective. In fact, in the pre-image of any such paraorthogonal 
polynomial is a measure that satisfies the hypotheses of the above theorems. Indeed, we 
have the following result, which is a consequence of the Bernstein-Szego Theorem (see [28 , 
Theorem 1.7.8]) 

Theorem 1.4. Let p be a probability measure with infinite support on the unit circle and 
suppose f3 e 90. Then 


®n(z-,P]fi) = 


g j(n—1)0 


dd\ 


<Pn-i(e i0 )(p^_ 1 (e ie ) 2t tJ 


where <p n -i(z) = <p n -i(z]fi). 


The following result will be relevant for our applications and is especially useful when 
combined with Theorem 11.41 


Proposition 1.5. Suppose w is as in Theorem 1.2 
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i) If w'{6) = f(e ld ) where f is analytic in a neighborhood of the unit circle, then the 
functions G n , D n , mid J n (defined in Theorem 1.2) can be analytically continued from 
{z : \z\ > 1} to {z : |z| > r} for some r < 1. 

ii) If w'{6 ) = f(e l6 ) where f is a rational function without poles on the unit circle, 
then the functions G n , D n , and J n (defined in Theorem 1.2) can be meromorphically 
continued from {z : \z\ > 1} to all of C as rational functions with no poles except 
possibly at zero and the poles of f inside D. 

Proof, i) We present the proof for G n \ the proof in the other two cases is even easier. We 
will make use of the fact that 


Therefore, we can write 

, 2 . ip* n _ l ( y e i6 )<p n _i(e ie )f {e i6 )_d6- 

2tt 


G n (z) = i 


f Tn-l(0<Pn-l(C)f(C)dC 


{z — e l6 )e ind 


(z - (K n+1 


2n' 


( 5 ) 


'o «/|C|—i 

If \z\ > 1, then our hypotheses allow us to move the contour of integration to the circle 
{C : ICI — t} for some r < 1 and hence we obtain an analytic continuation of G n (z) to the 
exterior of this circle. 

ii) Given the form of G n , D n , and J n , it suffices to show that if R(t) is a rational function 
of t without poles on the unit circle, then when \z\ > 1 

r 2lT R(e ie ) d.9 
J 0 z — e l6 27r 

is a rational function of z with no poles except possibly at 0 and the poles of R inside D. 
By using the partial fraction decomposition of R , this follows from the fact that if P is a 
polynomial and m is a non-negative integer, then 

P(e ie ) dd _ P( 0) 

J 0 z — e ld 27r z 

r 2n i in ( 1 


\z\ > 1, 


'0 


dO _ I z (- x )m 1 

(z — e l6 )(e l6 — x) m 27r | _ l _ 

v ' f z{z-x) m ' 


\z\ > 1, |x| > 1, 
\z\ > 1, |x| < 1. 


□ 


By combining the previous two results, we conclude that every paraorthogonal polynomial 
is a solution on all of C to a differential equation of the form (|2j) that has rational coefficients. 
This will be especially relevant for our applications. 

In the next section, we will discuss an application of polynomial solutions to linear second 
order differential equations. In Section [3] we will prove Theorems L2 and L3 The key idea 
will be to use the Szego recursion in several places to simplify various formulas. Finally, in 
Section [4] we present some detailed examples that highlight the applications discussed in the 
next section. 
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2 . Applications 

In this section, we will highlight some important consequences and applications of the 
main results of the previous section. 

2.1. Electrostatics. Orthogonal polynomials have proven to be useful tools when studying 
the equilibrium positions of electrons that repel each other via the two-dimensional Coulomb 
interaction. It was Stieltjes work that lead to the discovery that if n > 3 particles are confined 
to the interval [—1,1] interacting by means of a logarithmic (i.e. Coulomb) potential, then 
the unique energy minimizing configuration consists of one particle located at 1, one located 
at —1, and the other n — 2 particles located at the (distinct) zeros of a particular Jacobi 
polynomial (see [32], 03, 03] and see [36] for a proof). The key step in the proof of this 
fact is rewriting the equilibrium condition on a collection of points as a second order linear 
differential equation that is solved by the polynomial with zeros at precisely those points 
and then recognizing that a Jacobi polynomial with the appropriate choice of parameters 
is a solution to this differential equation (see [36] for details). Further developments in the 
electrostatic interpretation of zeros of orthogonal polynomials on the real line can be found 

m mmmmmm- 

Our main application of the results in Section [l] will be to demonstrate that the zeros of 
certain families of POPUC are the locations of points that are in electrostatic equilibrium. 
We will use the convention (as in m) that if a particle of charge q is located at a point 
a G C and a particle of charge p is located at a point b e C, then the force on the particle 
at b due to the particle at a is 2 pq/(b — a). 

We will consider the problem of creating an electric field that will keep identical charges at 
fixed points on the unit circle stationary. More precisely, suppose we are given a collection 
{aq,... ,x n } C <90. We will demonstrate a way to find a number m, a collection of points 
C C\{aq,... , aq}, and a set of real charges {qi\™ =l so that if a particle of charge 
+1 is placed placed at each aq (j = 1, ...,n) and a particle of charge g* is placed at a* 
(i = 1,... , m), then the total force on the particle at each aq is zero (j = 1 ,,n). In 
other words, if we have a collection of identically charged particles all lying on a concentric 
circle, we will demonstrate a way to construct an electric field that will keep these particles 
stationary. The points {a;}™ x and charges {g*}™i comprise what we call a set of electric 
field generators (in the language of [7j, the charged particles at {aq}” =1 would be called 
mobile charges and the charged particles at {a i \ 7 fL 1 would be called impurity charges ). It 
is important to keep in mind that the charged particles at the points {aq}” =1 interact with 
each other as well as with the electric field generators. 

With our motivation now clearly stated, we provide the following definitions: 
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Definition, i) Given a set of electric field generators {aj}™ 1 and {f/,}™ u we will say that 
a collection of points {aq,... ,x n } located on a smooth curve T is in T-normal electrostatic 
equilibrium if for each j — 1,... ,n, the force at aq is normal to T at Xj. 

ii) Given a set of electric fields generators and {qj}^, we will say that a collection 

of points {aq, ..., x n } C <90 is in total electrostatic equilibrium if for each j = 1 ,..., n, 


£ 

k =1 


1 


Xj Xfo 




Xj 


0. 


( 6 ) 


Notice that if Y is the unit circle, then the T-normal electrostatic equilibrium condition 
can be rewritten as 


Im 


/ 


Xj 


\ 


V-+ v 

t ^ O'* . _ o*, / ^ 


Qi 


— J Xj — Xh ' Xj — a, 
\ k=1 0 j=l 3 1 

\ 


= 0 , 


j = 1,2, 


(7) 


so it is clear that a collection of points on the unit circle that is in total electrostatic equi¬ 
librium is also in <9D-normal electrostatic equilibrium, but the converse is false. Indeed, it is 
well-known and easy to show that n particles of identical non-zero charge placed at the n th 
roots of unity and subject to no external force are in 9D-normal electrostatic equilibrium, 
but are not in total electrostatic equilibrium. 

The formula ([7]) is a restatement of the fact that if we place particles of charge +1 at each 
Xj (j = l,..., n) and a particle of charge g* at each a, (i — 1,..., m), then the condition ([T]) 
is satisfied if and only if the force exerted on the particle at Xj is normal to the unit circle 
at Xj for each j = 1 ,,n (see mm)- Similarly, the condition ([6]) is satisfied if and only 
if the force exerted on the particle at Xj is equal to zero for each j — 1,..., n (see [23]). 

As in the case of the interval, orthogonal polynomials are a useful tool when studying 
electrostatic equilibria on the circle. In [B], Forrester and Rogers use Jacobi polynomials to 
find a collection of 2 n points on the unit circle that is in <9B-normal electrostatic equilibrium 
when the electric field generators consist of a particle with charge p at 1 and a particle with 
charge q at —1. Their result assumes the added symmetry condition that each arc of the 
unit circle connecting 1 to —1 contains an equal number of points. 

Our main application is the following result: 


Theorem 2.1. Given any collection of n > 2 distinct points {aq,... ,x n } C <90, there exists 
a set of electric field generators so that the collection {aq,... ,x n } is in total electrostatic 
equilibrium. 


In fact, given any n distinct points {oq,...,x n } C 90, we can write down an explicit 
algorithm for finding the electric field generators. We proceed as follows: 

• Step 1: Define the measure p n on 90 by 


pn 


1 

n 



3 = 1 











and define f3 := (—l) n+1 YYj=i%j- 

• Step 2: Perform Gram-Schmidt orthogonalization on the linearly independent set 
{l,z ,... ,z n_1 } in L 2 (<9B,/i ra ) to arrive at the sequence of orthonormal polynomials 

+l(+ P"n)i • • • i V+i—1(+ /+)}• 

• Step 3: Define the probability measure 


dv n \ = 


d6 


|V3 n _i(e ie ;/i n )| 2 27r‘ 


• Step 4: Calculate the quantity h n (z ; f3] (3 ) for the measure v n in the domain {z : \z\ > 
1}. By Proposition [L5j it will be a rational function Si(z)/S 2 (z) for some polynomials 
Si and S 2 - 

• Step 5: Place a particle of charge —1/2 at each zero of Si, a particle of charge +1/2 
at each zero of S 2 , and a particle of charge 1(1 — n) at zero. 


The validity of the above algorithm will follow from the proof of Theorem 2.1 In order to 
prove Theorem 2.1 we need to translate the equilibrium problem into a differential equation 
so that we may apply the results of Section [TJ The appropriate differential equation in this 
setting is called a Lame equation, which we now discuss. 


2.2. The Lame Equation. A generalized Lame differential equation is a differential equa¬ 
tion of the form 


y'\z) 



y\z) 



— tV{z) = 0 , 

Wj) 


( 8 ) 


where S(z) is a polynomial of degree at most m — 2. Suppose we can find a polynomial P(z) 
of degree N with distinct zeros that solves ([8]) and also satisfies P(wy) ^ 0 for i — 1,..., m. 
Let be the zeros of P. It is easy to verify that 


so setting y = P and 2 = p 3 in Q shows 


N 

E 




Pk 


E 


ti /2 


= 0 


I =iPi 


Wi 




We recognize this as the condition for total electrostatic equilibrium with external field 
generated by a collection of charged particles located at where the particle at Wi 

carries charge U/2. We will see that the proof of Theorem [2T] follows from applying the above 
reasoning to the differential equation in Theorem |1.2| under the appropriate hypotheses. To 
this end, Proposition 1.5 will be essential. 

There is an extensive literature on the topic of generalized Lame differential equations and 
their relevance to electrostatics. We refer the reader to the references 0 EOl [2U [2211231 [251126] 






















for further information. We also refer the reader to [a El El HU] for further results concerning 
applications of polynomial solutions to second order differential equations. 


3. Proofs and Calculations 

Proof of Theorems M and \1.3[ Our proof begins very much in the same spirit as the proof 
of [16} Theorem 2.1]. We begin by writing 


ra-l j-2-K _ 

K&P) = ^Vk{z) / & n (e l9 ;/3)(p k (e ie )w(6) 

n JO 


k=0 

n —1 


de 

•a 


1 ‘ 2 -rr 


do 


^Lp k (z) I ie l9 & n (e %9 ] /3)e ie p k (e i9 )w(9) — 


k =0 
n—1 


/ 27r j _ jn 

— [$n(e zd ; P)]e ie ip k (e ie )w{d) — 

K = U 


We then integrate by parts to rewrite this as 


n 1 t> 27r j 

K&P) = -^2<Pk(z) J <$>n{e ie -p) — [e ie ip k (e ie ) 


k =0 
n —1 




f 2 IT 


/ $ n (e ie ]/3) e i9 (p k (e i9 )w'(d) - iw(8)e 2ie (p' k (e ie ) + e i9 (p k (e ie ) 


k =o 


dd 

2tt' 


where we used the continuity properties of w. For each m G N, let us define K m (z,t ) : = 
J2"j l =oTj( z )Tj(t) be the reproducing kernel for the measure /a and polynomials of degree 
at most m. By in Section 2.2], we can rewrite the above expression as 


in 

«h n (e^;/3)W„_ 1 (z,e^)e- i V(0)^ 

rZ7r _ _ jn 

+ (p n -i{z) J (e i9 tp n _i(e ie ) + e 2 V„-i(e*)) $n(e* e ; /3)w(d) — 

rlir in 

= i / ^ n ( e ie ^) K n _ 1 { z , e ie ) e - i 6 w \ 6 )— + n $ n - 1 { z){l - (9) 

Jo 

It follows from [37} Section 2.3] that (with ( = e l6 ) 


Kn-l(zX) 


l^n— lTn— 1 (.C) P<f) 

~z-C 


where 


ft 


C$n-l(Q 
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If we plug this into (9) and simplify (using also the fact that <p n -i( C) — C 1 ’Vn-i(C) when 
ICI = 1), we get 


& n (z;/3) = $ n _i(z) (n(l- pa n ^) + izK n - 1 


f 2 * ^ n {e ie -p)w'{d)^ n _ l ^ e ) dd 
l Q {z — e l6 )e ind 27r 




^ ^> n {e id ] p)w\6)^ n _ i(e* 0 ) dd 


J o (z — 27t 

= $ n -i(z)(n(l - j3a n - 1 ) + zG n (z) + j3D n (z)) - <h*_ 1 (^)(J n (^) - /3G n (z)), (10) 


One can also derive (10) from (J9| by using the Christoffel-Darboux formula (see m Section 
3]) to replace K n _\{z\ e l6 ). 


The relation (|10|) holds for all values of (3 G C. Let r G C be distinct from j3. Equation 

P®n(z-,r) -f$ n (z-,P) , v $n(z\P) -®n(z',T) 


Q easily implies 


$ 71 - 1 ( 2 ) = 


z{P - r ) 


K-i(z) = 




If we ping this into (10), we get the following system of ODE’s: 


rtf a\ * t \ /^M 1 ~P a n-i) +zG n (z) +(3D n (z)) t (J n (z) -pG n (z))\ 

= - W ~ F) -+-—§-J ( n > 

ml ■ M 1 “K-0 +zG n (z) + pD n (z)) | {J n (z) -pG n (z))\ 

n(z ’P J l z (p — t) f — B J 


rtf ^ ^7 \ p(n(l - TOn-i) + zG n (z) + rD n (z)) {J n (z) - rG n (z)) 

$nW T ) = <& n \Z\T) | - 7 ^-Z7-1" 


z{(3 - T ) 


T-(3 


( 12 ) 


_ $ ^ | V(w(l - ra B _i) + zG n (z) + tD u (z)) + (J n (^) - tG„(;z)) 


z{P - t) 


r-p 


This completes the proof of Theorem 1.3 We continue with the proof of Theorem 1.2 


Equation (11) yields the following formula for <3> n (^;r): 


z(P -T)& n (z;/3) + ® n {z-,P) [r{n( 1 - /3a n _i) + x;G n (^) + (3D n (z)) - z(J n (z) - j3G n {z))) 


P(n(l - (3a n -i) + zG n (z) + pD n (z)) - z(J n (z) - /3G n (z )) 


(13) 


We can use (13) to write 


K(z]t) = ( &^(z]P)z(P - f)h n (z;P;P)+ 

K(z;P) {(P - f)W[h n (z] /3; 0), z\ + h n (z;P;P)h n (z;r;p)) (14) 

+ $ n (z; P)W[h n (z; P\ P), h n (z ; r; /3)]^) h n [z\ P\ P)~ 2 
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If we substitute (14) into (12), then the resulting differential equation is 

( 4> “ (s; /J)z( ' 5 “ f)K(z ' /,) + ~ ft ft> -1 

+ h n (z;P;P)h n (z-,T-,P)) + $ n (z;/3)W[h n (z; 13; /3), h n (z-, r; /3)]f (15) 

_ h n (z;p;T)(z(fi-f)& n (z;p) + $ n (z; p)h n (z]T; p)) ^ Q ^ u 

7 / a a\ ®n\Zi P)h n {Z, T , TJ. 

hn\Z\P\P) 

This ODE holds for all r € <90 (even at P if \P\ = 1), so we can integrate both sides of this 
differential equation around 90 with respect to dr/(2ni). To do so, we will use the fact that 
f — 1/t and apply the Cauchy Integral Formula. Performing this lengthy calculation gives 
the ODE: 

-2 pW[h n (z;P;P),z\ 


0 2 Pz^(z; P) + z& n (z\ 0) ^ h n (z; P\ 0) 

+ P (n( 1 + (3a n _i) - j3D n (z ) - zJ n {z)/P) 
-zW[h n (z; P; P), h n (z ; ~P; P)\ 


- h n (z ; ~P;P) 


+ P) 


+ h n (z] P; P) [n + zG n (z)] 


h n (z;P;P) 

P ( (n{l - Pa n - 1 ) + £G n (£) + PD n (z )) (n - zJ n (z)/j3 ) 
T z{^J n {^Z) PC’ni.z'f) (uQ( n _] D n (z)) 


After some simplification, this can be rewritten as ([2]). □ 

In the proof of Theorem 1.2] the calculations up to equation (10) are very similar to those 
in the proof of pn Theorem 1.2]. However, our approach involving the system of first order 
equations given by 0 and ( p~2] ) differs from the approach used in [16] , which involves raising 
and lowering operators. In the next section we will see how these differing approaches result 
in different differential equations (even in the case P = a n _i), exemplified by the fact that 
the differential equation we derive depends on whether one chooses |z| < 1 or |z| > 1 when 
evaluating the required integrals. 


Proof of Theorem 2.1. First recall from [251 Theorem 2.2.13] that any collection of n distinct 
points on the unit circle is the zero set of a paraorthogonal polynomial on the unit circle. 
By Theorem 1.4] we may assume that the measure of orthogonality satisfies the hypotheses 


of Proposition |1.5[ i. 

If /i is a measure that satisfies the hypotheses of Proposition 1.5 i, then all of the coefficients 
in the differential equation in Theorem 1.2 are rational functions. Therefore, the ODE (|2]) 
can be written 


0 = <(*;£) + 


1 -n _ h' n (z] P\ P) 
z h n (z; P\ P) _ 


Kte ft + TrrKfeft, 


(16) 
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for some polynomials Qi(z) and Q 2 (z). Since h n (z ; /3; (3) is also a rational function, we can 
write it as S\(z)/S 2 (z) for some polynomials S\ and S 2 . With this notation, we can rewrite 
(16) as 


o = K(z;P) + 


1 ~ n _ S[(z) S' 2 {z ) 

2 Si(z) S 2 (z)_ 




If {kj t m}m =are the zeros of Sj (j = 1, 2), then we have 


0 = K(r,P) + 


deg(Si) 

1 — n v - 1 

z ^ z — k lr 

m =1 ’ 


deg(S 2 ) 


E TV 


m= 1 


z — k 


2,m 


Q 2 (z) 




(17) 


Once we show that $ n (z; j3) and Q 2 (z) do not share any common zeros, then our discussion 
in Section [272] implies the desired result. 


Notice that since p satisfies the hypotheses of Proposition |1.5f i, we can rewrite (16) as 

1 — n h' n (z;/9;j3y 


0 = $Z(z;/3) + 


z h n (z; /3; (3) 


®'n{ z 'i P) + 


Ri{z)h' n (z; P; P) 
hn(z;fcP) 


+ R 2 \Z) 


M z ;P), (18) 


where R\ and R 2 are rational functions without poles on the unit circle. Now, suppose for 
contradiction that Q 2 has a zero z 0 G <90 that also satisfies & n (z 0 ] (3) = 0. It follows that 
Q 1 /Q 2 das a pole at z 0 , which implies h' n (z ; f3; / 3)/h n {z ; f3\ (3) has a pole at z 0 (since Ri and 
R 2 have no poles on 50). Notice that h! n {z\ f3 ; f3)/h n {z ; /3; (3) must have a simple pole at zq 
and hence 

’ Ri(z)h' n (z] /?; P) 


1 ( O. o\ + R2\z) 

hn\Z) Pi P) 

has a removable singularity at zq. However, 

~ 1 ~n _ h' n (z)P)P) 
z h n (z; P‘, P) 


®n(z-iP) 




has a simple pole at zo (by the Gauss-Lucas Theorem), and hence the right-hand side of (16) 
is a rational function with a pole at z 0l which means it is not the zero function. This is our 
desired contradiction. □ 


We conclude this section by providing an explicit justification for the algorithm for finding 
the set of electric field generators outlined at the end of Section [2j From [2H, Theorem 
2.2.13], we know that with p n and P defined as in Step 1 of the algorithm, the polynomial 
z$ n _i(z] p n ) —PQ^-Az'i Pn) vanishes precisely at {ay,..., x n }. The Bernstein-Szego Theorem 
implies <F n _i(v; p n ) = $ n _i (z‘,u n ) {y n is defined in Step 3 of the algorithm). Therefore 

z<f> n _i(z;/x n ) - i3®n-i(z; Pn) = ®n(z',P', z'n)- 


Since v n satisfies the hypotheses of Proposition 1.5[ i, we can make the conclusion as in Step 
4 of the algorithm, and © shows that & n (z] /3; is n ) solves an ODE of the proper form to 
justify the placement of charges as in Step 5. 
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4. Examples 

In this section, we will consider applications of the results in Section [2] to specific proba¬ 
bility measures on the unit circle. 

4.1. Example: Lebesgue Polynomials. Let n be Lebesgue measure on the circle. In this 
case a n -i = 0 and w is constant so w' = 0. Let us also assume /3 = 1, so that ^.(z; (3) = z n —1 
and h n (z; <9; (3) = n. With this choice, the right-hand side of (J2]) becomes 

1 — n\ 

— H 

exactly as predicted. We see that particles of charge +1 located at the n th roots of unity are 
in total electrostatic equilibrium when the external field is generated by a charge of |(1 — n) 
located at the origin. 


n(n - 1 )z n ~ 2 + nz n_1 


4.2. Example: Degree One Bernstein-Szego Polynomials. In this case, let us write 


dfi(9) = 


l-ICI 2 dd 


ICI < 1. 


jl - (e ie \ 2 27t’ 

Some properties of this measure are given on |28], page 85], but there are some typos in the 
information there, which we will correct. If one sets £ = re llf and defines the Poisson kernel 
by 

i ji 

P r (x,y) := 


1 + r 2 — 2 r cos(a: — y) 

as on [211 page 27], then it holds that dy(9) = P r (9, — From this, it is easy to verify 
by direct computation that the Verblunsky coefficients satisfy «o = C an d a n = 0 for n > 1. 
The orthonormal polynomials {(p n (z; /i)} n >o f° r this measure satisfy 

We can explicitly calculate the ODE satisfied by & n (z;(3), though for notational conve¬ 
nience, we will specialize to the case ( = 1/2. The computations are lengthy, but each step 
can be handled using only simple contour integration and Fourier expansions. Indeed, if we 
assume \z\ > 1, then we calculate: 


GJz) = 


D n {z) = 


2(.~ 2 - 1 ) 


z(2z — 1) ’ ny j (2z — l) 2 ^ n_1 ’ 

With this knowledge of G n , D n , and J n , we can write 


J n {z) = 0. 


KM =^-D a +^-i)f; + ^-i) =; fjz) 


z n (2z — iy 


P2,n(zY 


where 


Pi,n(z) = z n (2z - l){n(2z - 1) + 2) + 2 zfi(z 2 - 1), P 2>n (z) = z n (2z - l) 2 , 
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If we write 


n +1 


Pl,n( z ) = Anz H(^ 

3 = 1 


then we can apply formula ( ]~Tt| ) and write 

2 


0 = & n (z;0) + 


n -\-1 


z — 1/2 


~[ Z ~ P3,r 


Pj: 








(19) 


where Qi.n and Q 2 >n are polynomials. Theorem 2.1 implies that if we place a particle with 
charge +1 at 1/2 and a particle with charge —1/2 at each point {pj, n }j=i, then n particles 
with charge +1 located at the zeros of <3> n (z;/3) will be in total electrostatic equilibrium. 
Figure 1 is a Mathematica plot illustrating this phenomenon when n = 22 and /3 = —1. 





♦ 



Figure 1. The squares show the location of the zeros of <3 >22 (2; —1) and the 
circles show the non-zero zeros of Pi^iz) when /3 = —1. The point {1/2} is 
marked by a diamond. 

Notice that the functions G n (z), D n (z), and J n [z ) are all holomorphic on D also. If we 
assume \z\ < 1 and n > 2, then we calculate 

1 2 z n ~ 2 (z 2 — If 

Gn(z) = D n (z) = 0 , J n (z) = ~ ^ j . 

With this knowledge of G n , D ni and J n , we can write 

_ 0z(z - 2 )(n(z - 2) - 2 z) - 2 z n (z 2 - 1) _ Pi,„(z) 

^-2)2 “'PynW’ 

where 


-Pi.n(^) = fiz(z — 2)(n(z — 2) — 2z) — 2 z n (z 2 — 1), 


P2,n(z) = z(z - 2) 2 , 
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If we write 


71+1 

P 1,71 (•2') 2^ J^J(^ Pj,n)i 

3 = 1 


then we can apply formula (JTt]) and write 


0 = K(z\P) + 


2 — n 


71+1 


+ 


z — 2 


E- 


3 = 1 


Pj,r 


®n( Z: >P) + 


Ql,r 


Q“2,r 




( 20 ) 


where Q liU and Q 2 , n are polynomials. Reasoning as above, we conclude that particles of 
charge +1 located at the zeros of <h n (z; /3) are in total electrostatic equilibrium when subject 
to the external held generated by a charge of +1 at 2, a charge of —1/2 at each pj >n , and a 
charge of |(2 — n) at the origin. Figure 2 is a Mathematica plot illustrating this phenomenon 
when n — 22 and j3 = —1. 


• ■ 


■ • 


■ • 


■ • 
■ • 


FIGURE 2. The squares show the location of the zeros of $22 ( 2 ; —1) and the 
circles show the non-zero zeros of Pi^iz) when f3 = —1. The point {2} is 
marked by a triangle and the origin is marked by a diamond. 


Since a charged particle at the origin exerts a force on a charged particle on the circle that 
is normal to the unit circle at that point, then if we remove the charged particle at zero from 
this example, the zeros of $„++) are still in <90-normal electrostatic equilibrium. 


4.3. Example: Sieved Degree One Bernstein-Szego Polynomials. In this case, we 
consider measures of the form 


dp{6 ) 


1-ICI 2 do 

1 - (e iMe \ 2 2vr’ 


ICI < l, MeN. 
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For this measure, the Verblunsky coefficients satisfy chm-i — C and a n — 0 for n ^ M — 1 
(se -Qb page 84]). Again, for the sake of clarity, we will specialise to the case of £ = 1/2 
and derive the second order ODE for $ n (^;/3). The calculations are very similar to those 


in Example 4.2, so we present fewer details here. The only additional tricks we use are the 
formulas 

„ M 


E 


x M — 1 


M 

E 

3 = 1 


£(£ M ~ x ) Mx \C, — x 1 / M e 2m i/ M £/ ’ 

M 1 


/M ^ ’ 


which are easily checked. When |z| > 1 and n > M, we have 

M _ , . 2M(> 2M -1) 


G„(*) = 


z(2z M - 1)’ 


D„U) = 


(2^ M - l) 2 z n ~ M ' 


Jn(z ) = 0 , 


h n (z;/3; /3) = (3 ( n + 


2M 


+ 


2 M(3(z 2M - 1) 


2z M -l {2z M -l) 2 z n ~ M / 

If we perform an analysis similar to that of the previous example, then we deduce the 
existence of a polynomial P\ >n {z) given by 

n+M 

Pl,n( Z ) '■= CzM II(- 


for some C G C so that & n (z] /?) satisfies 

M 


3 = 1 


n+M 


KP, + E — 


0 = 1 




1-M 1 

H-/ - 

- jzf * - Pp 


V2,nbJ 


where Qi,n and Q 2 ,n are polynomials. Therefore, particles located at the zeros of 4> n (z;/3) 
- each carrying charge +1 - are in total electrostatic equilibrium when the external field is 
created by a particle of charge +1 at each of {2 a particle of charge —1/2 at 
each of {pj, n }j=, and a particle of charge 1(1 — M ) at the origin. If we remove the charge 
at the origin, then the zeros of < h n (^; {3) are in <9D-normal electrostatic equilibrium. 


4.4. Example: Single Non-Trivial Moment. In this case, we consider the measure given 
by 

J t q\ f, ,a\\ de |1 — e u '\ 2 d6 
dp{6) = (1 - cos(0))— =---—. 

For this measure, the Verblunsky coefficients satisfy a n = — (n + 2) _1 (see (28], page 86]). 
The monic orthogonal polynomials are given by 


®Jz) = 


n +l z 

3=0 


Eo +a 


1 As in the previous example, we use the table on page 85 in [28] with P r (Q,ip) replaced by P r (9 , —ip). 
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From this, it follows that 


$n(z) “I) = 


„ra+l 


which has zeros located at the (n + 1+ roots of unity, except one. This means that particles 
of charge +1 located at the zeros of <3> n (z; —1) are in ciD-normal electrostatic equilibrium 
when the external field is generated by a single particle of charge +1 located at 1 (see also 
[6] or |H Theorem 3] with p — 1 and q = 0). 

With the above formulas, and the fact that K n = we calculate (for \z\ > 1) 

1 


G n {z) — - / | -i n 

n(n + 1 )z 


" k 2 

E L 


n — n 


D n {z) = 

Jn(z') 


-1 


n(n + 1) 

-1 

zn(n + 1) 


fc=i 


n 


E 

fc =i 


n 2 + 2 n — 2 kn — k 2 


n 


yU +1 


n—1 

E 

,k =0 


(n — k)" 


With this knowledge of G n , D n , and ./„, we can write down h n (z]/3] (3) as a ratio of two 
polynomials, but the expression is extremely lengthy. For notational convenience, we consider 
only the case j3 — —1: 

n{nz n+2 — {n + 2)z n+1 + z + 1) Pi,n(z) 


hn(yZ] 1; (n + l)z n+1 (z — 1) 

where 

Pi,n(z) = ~n(nz n+ 2 - (n + 2 )z n+1 + z + l), 

If we write 

Pi,n(z) = -n 2 {z - 1) JJ( 


' P2,n{zV 

Pi,n{z) = (n + 1 )z n+l {z - 1). 


71+1 


Pjr 


3 = 1 

then the differential equation of Theorem L2 becomes 

n+1 


0 -<+-!)+ ( --Y.VZ [ 


3 = 1 






where Qi,n and Q 2 , n are polynomials. Reasoning as above, we conclude that particles of 
charge +1 located at the zeros of < F n (^;—1) are in total electrostatic equilibrium when the 
external held is generated by particles of charge —1/2 at each p j n and a particle of charge 
+1 at the origin. Figure 3 is a Mathematica plot illustrating this phenomenon when n = 14. 

For the sake of comparison with the results in [I6j, let us compute the differential equation 


from Theorem 1.2 when setting j3 = a n _i so that Q n (z]f3) = < F n (^). This example was 

ere that 

+ &n(z)^z -r = 0. 


considered in PH Example 1] and it was shown there that the polynomial < h n (^) satisfies 

2 n 


+<(*) (-3 - 


zil-z) 
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FIGURE 3. The squares show the location of the zeros of $i 4 (z; —1) and the 
circles show the zeros of -Pi, 14 (F) except z — 1. The diamond is located at the 
origin. 

In our calculations, we find (using Mathematica to simplify the expressions) 

— 1 —1 \ n(z n+3 (n + l ) 2 — z n+2 (2n 2 + 6 n + 3) + z n+1 (n + 2 ) 2 — z — 1) 

z n+1 (n + l) 3 (z — l ) 3 

(. z n+2 (n + 1 ) — z n+1 (n + 2 ) + 1 )(z(n + 1 ) + {n + 2 )) 


+ 


h n Z] , 1 , , , 

One can then calculate 

^ l n { Z i n+1’ n+l) _ 1 . 

(z; z 1 — z' z(z n+3 (n + l ) 2 — z n+2 (2n 2 + 6 n + 3) + z n+1 (n + 2 ) 2 — z — 1) 

and observe that the differential equation we derive is different than the one found in 
As mentioned in Section [lj this distinction is not unexpected. 

If we define G n , J n , and D n for z e O, then for n > 2 we have 

”n—1 

y: (^ fc (n 2 + n — (k + l) 2 )) 


= 


-1 


n(n + 1) 


7 ' 77 

n(n + 1) 


./c =0 

n—2 


Jyi.(z) 


-1 


n(n + 1) 


k =0 
2 


53P*(fc-3n+l) 2 ) 

0 

n—1 

z n + (**((* + n + l ) 2 - 2 n - 2 n 2 )) 


fc =0 


With this knowledge of G n , D n , and J n , we can write down h n {z\ (3; (3) as a ratio of two 
polynomials, but the expression is extremely lengthy. For notational convenience, we consider 
only the case j3 — —1: 

n(z n+2 + z n+1 — (n + 2)z + n) P ljn (z) 


hn^Zi 1 ) 1 ) 


(n + l)(z - 1) 


" P2,n(*Y 
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where 


1,71 

If we write 


Pi.Jz) = n(z n+2 + z n+l - (n + 2 )z + n), P 2 ,n(^) = (n + l)(z - 1), 


n+l 


Pi,n{z) = n{z - 1) - p j>n ), 

3 =1 


then the differential equation of Theorem L2 becomes 

0 = tlfe -1) + ( i-2 - £ -i_) t ; (2 ; -1) + -l), 


J=1 


Pj* 


Q2,n (^) 


where Qi,n and Q 2 ,n are polynomials. Reasoning as above, we conclude that particles of 
charge +1 located at the zeros of 1) are in total electrostatic equilibrium when the 

external held is generated by particles of charge —1/2 at each pj jU and a particle of charge 
|(1 — n) at the origin. Figure 4 is a Mathematica plot illustrating this phenomenon when 
n = 14. 


FIGURE 4. The squares show the location of the zeros of $14 (z; —1) and the 
circles show the zeros of Pi 14: (z) except z = 1. The diamond is located at the 
origin. 
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